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Finite volume corrections to higher moments are important observable quantities. They make 
possible to differentiate between different statistical ensembles even in the thermodynamic limit. It 
is shown that this property is a universal one. The classical grand canonical distribution is compared 
to the canonical distribution in the rigorous procedure of approaching the thermodynamic limit. 
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Fluctuations and correlations measured in heavy ion collision processes give better insight into 
dynamical and kinematical properties of the dense hadronic medium created in ultrarelativistic 



^ , heavy ion collisions. Particle production yields are astonishingly well reproduced by thermal mod- 

els, based on the assumption of noninteracting gas of hadronic resonances 0]. Systems under 
considerations are in fact so close to the thermodynamic limit that final volume effects can be 
neglected — at least when productions yields are considered. 

The aim of the paper is to show that finite volume effects become more and more important 
^ . when higher moments, e.g. correlations and fluctuations are considered. The basic physical char- 

^\ ' acterization of the system described by means of the thermal model are underlying probability 

OO , densities that given physical observables of the system have specified values. The only way to 

reproduce those probability distribution is by means of higher and higher probability moments. 
Those moments are in fact the only quantities which are phenomenologically available and can be 
used for the verification of theoretical predictions. Finite volume effects are also important for the 
f-*) I lattice QCD calculations. 

Particle yields in heavy ion collision are the first moments, so they lead to rather crude com- 
parisons with the model. Fluctuations and correlations are second moments so they allow for the 
, better understanding of physical processes in the thermal equilibrium. 

A preliminary analysis of the increasing volume effects was given in 0, Q • It has been rigorously 
(— I I shown an influence of 0{1/V) terms for a new class physical observables — semi-intensive quan- 

tities 0. Those results completely explained also ambiguities noted in [3], related to "spurious 
non-equivalence" of different statistical ensembles used in the description of heavy ion collision 
processes. 

' This paper is devoted to a further analysis ofO{l/V) terms. It is shown that those terms are not 

5^ , specific for systems with subsidiary internal symmetries but appear also in the simplest "classical" 

problems of statistical physics. 

II. CHOICE OF VARIABLES 

In the thermodynamical limit the relevant probabilities distributions are those related to densi- 
ties. These distributions are expressed by moments calculated for densities — not for particles. In 
the practice, however, we measure particles — not densities as we do not know related volumes. 
Fortunately, volumes can be omitted by taking corresponding ratios. 
Let us consider e.g. the density variance An^. This can be written as 
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By taking the relative variance 

An2 _ (jV2) - (jV)2 

volume-dependence vanishes. 

A. Semi-intensive variables 

A special care should be taken for calculations of ratios of particles moments. Although moments 
are extensive variables their ratios can be finite in the thermodynamic limit. These ratios are 
examples of semi-intensive variables. They are finite in the thermodynamic limit but those limits 
depend on volume terms in density probability distributions. One can say that semi-intensive 
variables "keep memory" where the thermodynamic limit is realized from. 

Let consider as an example the scaled particle variance 

{N^)^{N)^ _ {n')-{nf 
{N) in) ■ 

The term 

¥) ■ 

tends to zero in the thermodynamic limit as 0{V^^). So a behavior of the scaled particle variance 
depends on the 0{V~'^) term in the scaled density variance. A more detailed analysis of semi- 
intensive variables is given in Q. 

To clarify this approach let us consider a well known classical problem of Poisson distribution 
but taken in the thermodynamic limit. 



III. GRAND CANONICAL AND CANONICAL ENSEMBLES 



A. Poisson distribution in the thermodynamic limit 



Let us consider the grand canonical ensemble of noninteracting gas. A corresponding statistical 
operator is 

D = — ^ (1) 

This leads to the partition function 

ZiV,T,j)^ e^^'' . (2) 
where z is one-particle partition function 

^(T, V) = / d^p e-'^^ff) = VzoiT) , (3) 



(27r)3 

A 7 parameter {— (3^) is such to provide the given value of the average particle number {N) = 
V{n). This means that 

e- = M. (4) 



Particle moments can be written as 
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The parameter 7 is taken in final formulae as a function 7((n), zq) from Eq (|4l). 
The resulting probability distribution to obtain N particles under condition that the average 
number of particles is (N) is equal to Poisson distribution 



We introduce corresponding probability distribution V for the particle number density n = N/V 



V^r.){n;V)^VPvin){Vn) = V 



T{Vn+l) 

For large Vn we are using an asymptotic form of Gamma function 

r{Vn + 1) - v^(T/n)^"-i/2 ^-Vn L ^ ^_ ^ 0(y-2)l 

12V n J 



This gives 



1 



12Vn 



0{V' 



(6) 



(7) 



This expression in singular in the V 00 limit. To estimate a large volume behavior of the 
probability distribution ([6]) one should take into account a generalized function limit. So we are 
going to calculate an expression 



(G)y = J dnG{n)Vi^n){n-V), 

where T't^n){'n] V) is replaced by the asymptotic form from Eq ([7|). In the next to leading order 
in 1/y one should calculate 



1 



(8) 



n-K J 12^/2^7 n3/2 

where 

S{n) — n ln(n) — n In n + n — (n) . 
An asymptotic expansion of the function (G) y is given by the classical Watson-Laplace theorem 
Theorem 1 Let I — [a, h] he the finite interval such that 

1. maxS'(a::) is reached in the single point x — xq, a < xq < b. 

2. f{x),S{x)eC{I). 

3. f{x),S{x) G C°° in the vicinity ofxg, and S (xq) 7^ 0. 
Then, for A — > 00, X d S^, there is an asymptotic expansion 



F[X] 



Ck 



k=G 

r(/c + 1/2) / d 

(2fc)! [jh: 



2k 



S{x^) ^ S{x )\-''-''^ 
[x - XoY 



(9a) 
(9b) 



is here a segment \ argz| 5j ^ — e < -I in the complex z-plane. 
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To obtain 0{1/V) formula the first term in ^ should be calculated till the next to leading order 
term in 1/V . For the second term it is enough to perform calculations in the leading order only. 
The first term gives the contribution 



and the second term gives 



So we have eventually 



1 / dn^ e^^(") = T^G((n)) , (10b) 
12\/27r7 n-^l^ l2{n)V 



{G)v = G{{n)) + i^G ((n)) + 0{V~^) , (11) 



for any function G. 

This gives us the exact expression for the density distribution Q in the large volume limit 

P(„) (n; V) ~ 5{n (n)) + M 5- [n - («)) + 0{V-^) . (12) 
We are now able to obtain arbitrary density moments up to 0{V~'^) terms. 

(n^V = Jdnn'^V^n){n;V) = {n)'^+'^^^^{n)'^~'+0{V-'). (13) 
We have for the second moment (intensive variable!) 



This means 

An^^M^O. (14) 

as expected in the thermodynamic limit. 

The particle number and its density are fixed in the canonical ensemble so corresponding vari- 
ances are always equal to zero. The result (jl4l) can be seen as an example of the equivalence of 
the canonical and grand canonical distribution in the thermodynamic limit. This equivalence is 
clearly visible from the Eq where the delta function in the first term can be considered as the 
particle number density distribution in the canonical ensemble. 

A more involved situation appears for particle number moments (extensive variable!). Eq p3p 
translated to the particle number gives 

(iV*^) + 1+0(^^-2), (15) 

One gets for the scaled variance (semi-intensive variable!) 

AA^2 

what should be compared with zero obtained for the canonical distribution. 

The mechanism for such a seemingly unexpected behavior is quite obvious. The grand canonical 
and the canonical density probability distributions tend to the same thermodynamic limit. There 
are different however for any finite volume. Semi-intensive variables depend on coefficients at those 
finite volume terms so they are different also in the thermodynamic limit. 
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B. Energy distribution 



It is interesting to perform similar calculation for the energy distribution in both ensembles. 
Energy moments and an average energy density can be written as 



{E') - (-1)' 



1 d'^Z 



One gets from Eq (fTT]) 

{E'') =V''{e)'' + V' 



2 ^ ' zo df3^ 



k-2\ 



The grand canonical energy density distribution follows 



P(6|(n>,(6»=<5(6-(6» + i^7^<=^ 



[n) 



5 (e-(6))+0(F-2), 



•^GG jg giygjj ]^epg as 



GC 



1 d^ZQ 



zq dp"^ 



/3=/3((e)/(n)) 



For the canonical distribution a corresponding statistical operator is 



D 



This leads to the partition function 



gVn log z 



m 



Internal energy moments are given by Eq (I17|) . In particular 

n dzQ 
^ ~An- 

Zq dp 

For the energy moments one gets now 
A corresponding probability distribution is 



Pie\n,{e))=S{e-{e)) + —n^ 



where TiP is given here as 



d_ / 1 dzp 



(17) 
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